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ABSTRACT 

The  research  summarized  concerns  several  aspects  of  the  propagation  and  scattering 
of  acoustical  and  optical  waves.  The  topics  discussed  fall  under  the  following  four 
categories: 

A.  A  generalization  of  the  geometrical  theory  of  diffraction  which  escribes 
surface-guided  wave  contributions  to  the  scattering  of  sound  from  smooth 
elastic  objects  in  water;  the  effect  of  curvature  on  Lamb  wave  velocities  on 
thick  shells;  a  product  expansion  of  the  S  function  for  elastic  spheres. 

B .  Acoustical  and  optical  caustics  and  associated  wavefields:  unfolding  axial 
caustics;  observations  of  lips  caustics  in  light  backscattered  from  oblate  drops; 
and  observations  of  acoustical  and  optical  transverse  cusps  produced  by 
reflection. 

C .  Physical  optics  of  bubbles  in  water:  observations  of  Brewster  angle  light 
scattering;  observations  of  backscattering  from  rising  bubbles;  asymptotic 
series  for  critical  angle  scattering  useful  for  determining  bubble  sizes. 

D.  Acoustic  phase  conjugation  and  aspects  of  wavefront  reversal:  a  theory  for 
wavefront  reversal  via  three-wave  mixing  in  a  layer  of  enhanced  nonlinearity  as 
caused,  for  example,  by  a  layer  of  microbubbles  in  water,  the  present  status  of 
experiments  on  three-wave  mixing;  and  the  possibility  of  phase  conjugation  via 
reflection  from  curved  vibrating  surfaces. 
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I.  SCATTERING  OF  SOUND  FROM  ELASTIC  OBJECTS  IN  WATER 

A.  Project  to  Extend  GTD  to  Smooth  Elastic  Objects 

The  extension  of  the  geometrical  theory  of  diffraction  (GTD)  to  backscattering  from 
elastic  objects  in  water  should  allow  complicated  scattering  problems  to  be  partitioned  into 
geometry  and  the  determination  of  local  guided-wave/acoustic-field  interactions  from 
continuum  mechanics.  It  gives  a  simple  and  quantitative  understanding  of  the  scattering 
process  which  could  be  useful  both  for  inverse  problems  and  for  predicting  how  changes  in 
an  object  will  affect  the  scattering.  Heuristic  support  for  the  basic  ideas  were  put  forth  by 
Soviet  researchers.1  In  their  approach,  however,  the  coefficient  which  describes  the 
coupling  of  surface-guided  waves  with  the  acoustic  field  was  determined  by  a  fitting 
procedure.  Independently,  and  at  about  the  same  time,  a  rigorous  analysis  of 
backscattering  from  elastic  spheres  was  carried  out  by  Williams  and  Mars  ton2-1  based  on 
the  Watson  transformation.  That  analysis  gave  an  exact  expression  for  the  required 
complex  coupling  coefficient  G/  for  the  specific  case  of  a  solid  elastic  sphere.  In  that 
expression,  however,  the  dependence  of  G /  on  physically  relevant  parameters  was 
obscured  by  the  complexity  of  the  expression. 

A  desirable  form  for  G /  is  one  in  which  it  manifests  its  dependence  on  the 
properties  of  the  surface  wave  on  the  object.  These  properties  are  the  local  phase  velocity 
ci  and  the  radiation  damping  parameter  P/.  Such  relationships  were  obtained  for  IG/I  in 
the  preceding  contract  year  with  the  following  result  for  the  sphere  and  circular  cylinder 
cases,  respectively4'5 

lG*pl  -  8:tPzc/cr  iGfl  =  87tP/(7tx)1/2,  (1,2) 

where  x  =  ka,  a  is  the  radius  of  the  object  and  k  =  2%fk  where  X  and  c  are  the 
wavelength  and  sound  speed  in  water.  These  expressions  were  confirmed  with  the  exact 
result  for  solid  spheres  and  Veksler's  fitted  result  for  Lamb  waves  c n  a  hollow’  cylinder.4*5 
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B.  Approach  and  Accomplishments  on  Extending  GTD 

One  unresolved  issue  was  the  phase  of  the  complex  coefficient  G /  for  the  various 
cases  of  interest.  These  were  calculated  for  certain  special  cases,  as  described  below,  by  a 
tedious  extension  of  the  method  used  in  the  derivation  of  Eqs.  (1)  and  (2).  The  approach  is 
based  on  the  observation  that  the  Fabry-Perot/GTD  representation  of  the  form  function 
should  reduce  to  the  result  of  the  Breit-Wigner  form  of  resonance  scattering  theory  (RST) 
with  an  appropriate  background  for  the  special  case  of  small  (}/  where  the  Breit-Wigner 
form  should  be  applicable.  While  the  resulting  approximations  for  IG/I,  Eqs.  (1)  and  (2), 
do  not  depend  on  the  choice  of  the  background,  the  resulting  approximations  for  the  phase 
do.  The  results6  are  summarized  below  where  the  superscripts  r  and  s  designate  the 
rigid  and  soft  background  choices,  respectively 

arg(G*p)r  -  0  ,  arg(G*p)s  -  K,  (3,4) 

arg(G^y)s  *  n/4,  arg^)*  «  5w/4 .  (5,6) 

where  sp  denotes  sphere  and  cy  denotes  circular  cylinder.  Equations  (3)  and  (5)  should 
be  applicable  to  solid  scatterers  which  are  sufficiently  dense  that  RST  with  a  rigid 
background  choice  is  a  good  approximation  to  the  exact  scattering.  Equations  (4)  and  (6) 
may  only  be  applicable  to  very  thin  shells  for  which  a  soft  background  causes  RST  to  be  a 
useful  approximation. 

The  correctness  of  the  above  analysis  was  confirmed  in  one  case  by  comparison 

with  exact  results  from  a  partial- wave  series  (PWS),  Fig.  1.  The  GTD-based  synthesis 

shown  uses  the  approximation  f  *  f$  +  f/=R  +  f/=WGi  where 
-2 (K-e,)p,  in, 

-G.  e  C  y  Pf^I  "  P^  \  -i2ka 

f,  =  - - - - - ,  G,  -  8rtp  c/c,,  f$  -  (JLL - )  e  2k  (7a,b,8) 

[l  +  j  exp(-2npi  +  i27cxc/c/)J  Pecl  +  Pc 

where  for  the  present  problem  of  spheres  j  =  +1.  Equation  (8)  is  a  simple  approximation  of 
the  specular  contribution  for  a  solid  sphere  of  density  pg,  longitudinal  sound  speed  cl  in 
water  of  density  p.  The  agreement  between  the  exact  and  GTD  synthesis  of  Ifl  in  Fig.  1  is 
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Fig.  1.  The  solid  curve  is  the  exact  If!  from  the  partial  wave  series  for  backscattering  from  a  tungsten  carbide  sphere  in  water. 
The  dashed  curve  is  the  simple  synthesis  based  on  Eqs.  (7)  and  (8)  as  described  in  the  text  The  computations  were  done  at  the 
Naval  Coastal  Systems  Center  in  collaboration  with  Kevin  L.  Williams. 
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satisfactory  considering  only  the  Rayleigh  wave  and  slowest  whispering  gallery  wave 
contributions  were  included  in  the  sum  over  /. 

In  a  related  accomplishment,  an  approximation  was  derived  in  a  heuristic  way  for 
the  contribution  of  a  surface  wave  to  scattering  from  a  circular  cylinder  not  restricted  to  the 
backward  direction.7  Furthermore,  it  was  possible  to  use  GTD  methods  to  relate  arg(G/sP) 
to  argCG/cy)  directly,  without  reference  to  RST.  It  was  assumed  that  ka  »  1  and  p/ 
and  c/  are  the  same  for  the  sphere  and  for  the  cylinder  so  that  the  targets  are  either  both 
solid  or  both  shells.  The  result  is  argCG/sP)  =  argufy)  -  (re/4)  which  agrees  with  Eqs. 

(3)  -  (6).  The  method  of  comparing  cylinder  and  sphere  scattering  is  like  that  of  Sec.  VI  of 
Ref.  5  except  that  now  the  phase  was  found. 

The  surf  ace -elastic-wave/GTD  approach  to  modeling  Lamb  wave  contributions  to 
backscattering  from  a  hollow  shell  was  confirmed  in  our  earlier  work  by  comparison  with 
data.4*8  Though  there  was  agreement  between  theory  and  data,  without  the  use  of  any 
adjustable  parameters,  it  was  necessary  to  use  the  Watson  transform  methodology  to  find 
the  Lamb  wave  parameters  c /  and  P/.  It  would  be  desirable  to  demonstrate  that  these  can 
be  found  from  the  local  continuum  mechanics  of  a  water  loaded  curved  plate. 

A  simple  approximation  was  obtained  which  is  an  advancement  toward  the 
aforementioned  goal.  In  this  approximation  the  effect  of  curvature  on  the  phase  velocity 
along  the  outer  surface  of  a  curved  plate  was  estimated.9  The  result  was  compared  with  the 
results  of  the  Watson  transform  for  the  case  of  a  spherical  stainless-steel  shell  having  a  ratio 
of  inner  to  outer  radius  b/a  of  0.838  and  a  =  19.05  mm.  The  shell  thickness  is  h  =  a  -  b 
=  3.1  mm.  In  the  discussion  which  follows  l  =  &o  designates  the  antisymmetric  (i.e. 
flexural)  Lamb  wave  while  /  =  sQ  is  the  symmetric  (i.e.  extensional)  Lamb  wave.  First 
consider  the  exact  results  from  the  Watson  transform  for  c//c  for  the  sphere  in  water.  The 
set  of  curves  which  begin  at  ka  =  20  in  Fig.  2  give  the  predicted  c//c  with  the  solid  curve 
being  l  =  ao  and  the  dashed  curve  being  l  =  s0.  The  curves  which  begin  at  ka  =  0  and  lie 
below  the  set  of  exact  curves  (described  above)  are  given  by  taking  c /  =  (c/)p  where  (c/)p 


o 

\  2 


’r  // 


Fig.  2.  The  set  of  curves  which  begin  at  ka  =  20  give  the  exact  c i/c  according  to  the 
Watson  transform  for  Lamb  waves  on  a  shell  in  water  where  c  is  the  speed  of  sound  in 
water.  The  lower  of  the  solid  and  short-dashed  curves  use  Lamb's  result  (c/)p  for  a  flat 
plate.  The  other  set  (which  are  in  much  better  agreement  with  the  exact  theory)  incorporate 
the  curvature  connection  given  by  Eq.  (9). 
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is  the  velocity  of  the  respective  wave  on  a  flat  plate  of  thickness  h  in  vacuum  as  given  by 
Lamb's  original  analysis.  For  this  set,  the  solid  curve  has  /  =  ao  while  the  dashed  curve 
has  /  =  Sq.  The  curves  with  the  alternating  long-short  dashes  and  the  long  dashes  are  a 
much  improved  approximation  in  which 

c,  =  (c;)p/[l  -  (h/2a)].  (9) 

The  denominator  in  Eq.  (9)  is  a  novel  method  of  accounting  for  the  curvature  of  the  shell. 
The  velocity  of  the  wave  at  the  midpoint  of  the  shell  material  [where  the  radius  is  a  -  (h/2)] 
is  taken  to  be  (c/)p;  this  velocity  is  projected  to  the  outer  surface  (of  radius  a)  to  give  c /. 
The  general  agreement  with  the  exact  result  shows  that  c /  could  have  been  obtained  from 
local  surface  properties  had  we  not  derived  the  exact  partial-wave  series. 

In  related  work,  a  GTD  based  phenomenological  model  was  examined  for  the  broad 
mid-frequency  enhancement  of  backscattering  from  thin  spherical  shells.  A  graduate 
Student,  Cleon  Dean,  participated  in  the  computations. 


C.  Product  Expansion  of  the  S  Function  for  Scattering  from  Elastic 
Spheres  Having  Multiple  Resonances 

In  related  work,  the  implicit  assumptions  of  formal  resonance  scattering  theory10 
have  been  corrected  and  clarified.  Associated  with  the  scattering  phase  shift  5n  of  the  nth 
partial  wave  for  a  sphere  of  radius  a  is  the  function  Sn(x)  =  exp[2i8n]  where  x  =  ka.  A 
theorem  from  classical  scattering  theory11  was  applied  by  Marston  to  obtain  the  following 


product  expansion  which  appears  to  be  important  and  novel  for  acoustics 


(x  i  -  x)(x  ,  +  x)  (iL  .  -  x) 

S  s  ± e~2  n  — ^ — -n7,nj 

n  ,  *  .  :  fit  ■+■ 


'  (X*  +  x)(x  ,  -  X)  J  (iLnj  +  X> 


where  the  xnt  =  Xn/  -  i(Tn//2)  lie  in  the  fourth  quadrant  and  Lnj  >  0.  Unlike  a  sum 
expansion  (stated  in  some  acoustic  RST  literature10)  Sn  remains  manifestly  unitary  even 


11 


for  multiple  resonances  /.  This  Sn  may  be  used  to  split  off  elastic  contributions  to  the 
partial-wave  form  function  fn  for  backscattering  in  analogy  with  RST.  In  the  case  of  only 
two  resonances  Qabeled  /  =  1  and  2) 


~f<b)  + 


nl 


+  fn2  + 


^int) 


(11) 


where  fn®  is  a  background  term  and  it  is  assumed  that  x  +  Xn/ »  r„/.  The  fn/  have  a 
Breit-Wigner  form 

fn/  =  exp(2i^)[(2n+l)/x](-l)n  T /[Xn/  -  x  -  (i/2)r/1  (12) 


and  is  the  phase  shift  associated  with  fn^X  The  interaction  term  vanishes  as 
Fn//lx  -  Xn/I  — »  0  for  /  =  1  or  2.  The  implicit  assumptions  of  RST  are  that  Sn  has  the 
product  form,  Eq.  (10),  and  not  the  sum  expansion  usually  written.10  Furthermore,  the 
interaction  term  ^(*"0  is  neglected  when  fn  is  written  as  a  sum  of  Breit-Wigner  terms. 

The  product  expansion  over  index  j  in  Eq.  (10),  manifestly  corresponds  to  poles 
(or  generalized  resonances)  for  which  there  is  no  restoring  force  but  only  damping.  (These 
poles  lie  on  the  imaginary  axis  at  x  =  -  iLnj.)  The  physical  significance  of  such  a  pole  is 
illustrated  by  considering  the  n  =  l  or  dipole  partial- wa*  x  The  Lij  pole(s)  account  for 
the  translational  motion  of  a  movable  elastic  sphere. 


II.  ACOUSTICAL  AND  OPTICAL  CAUSTICS  AND  ASSOCIATED 
WAVEFIELDS  (DIFFRACTION  CATASTROPHES) 

A.  Project  to  Explore  Novel  Acoustical  and  Optical  Diffraction 
Catastrophes 

In  previous  research  we  described  examples  of  acoustical  and  optical  diffraction 

i 

catastrophes  in  certain  reflection  and  scattering  problems.4  Since  such  catastrophes 

decorate  the  caustics  of  geometrical  optics,  the  amplitudes  of  the  wavefields  tend  to  be  large 

when  the  wavelength  is  short.  Three  recent  accomplishments  are  described  below. 

I 
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B .  Unfolding  Acoustical  Axial  Caustics  of  Glory  Scattering  with 
Harmonic  Angular  Perturbations  of  Toroidal  Wavefronts 

This  theory  and  an  analysis  of  the  associated  acoustic  wavefields  were  a  part  of  the 
recently  completed  Ph.D.  dissertation  project  of  W.  P.  Amott.12'13  Attention  was 
restricted  to  the  simplest  class  of  harmonic  angular  perturbation  f(y)  of  the  toroidal 
wavefront  W(s): 

W(s)  =  -i!^L,  f(¥)  =  4  (1+ cospy),  (13,14) 

2a  2 

where  s  and  \\r  are  radial  and  angular  polar  coordinates,  the  positive  constants  b  and  a 
each  have  dimensions  of  length  and  specify  the  focal  circle  radius  and  position  for  the 
unperturbed  wavefront  while  8  is  the  strength  of  the  perturbation  and  p  =  2, 3, 4  ... 
describes  the  harmonic  frequency  of  the  perturbation.  The  outgoing  wavefronts  of  interest 
are  of  the  form 

W(s,\|0  =  W(s)  -  f(\jr),  (15) 

which  is  relevant  to  several  scattering  or  reflection  problems  for  nonspherical  objects.12*13 
(One  example  is  Lamb  wave  contributions  to  backscattering  from  slightly  oblate  or  prolate 
shells  for  which  case  p  =  2.)  Let  the  z  axis  be  the  extension  of  the  axis  along  s  =  0 

where  W  =  0  corresponds  to  the  plane  z  =  0.  It  is  convenient  to  restrict  attention  to  weak 
perturbations  for  which  I  pip2  <  0.025  where  {3  =  8a/2b.  The  propagation  of  the 

wavefront  given  by  (15)  yields  a  caustic  surface  with  Cartesian  coordinates  (x,y,z)  which 
are  well  approximated  by 

x  =  (?..+  ^.pPb  {(p- l)cos[(p  +  l)y]  +  (p+  l)cos[(p- l)\jr]}  ,  (16) 

2a 

y  =  ■”.*  -pPb  { (p  -  1)  sin[(p  +  1)\|/]  -  (p  +  I)  sin[(p  -  lty]  }  .  (17) 

2a 

This  is  a  novel  result.  The  intersection  of  the  plane  of  fixed  z  with  the  caustic  surface 
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traces  out  the  caustic  curves  shown  in  Fig.  3  where  p  =  2, 4, 6,  3,  and  5  in  (a)  -  (e), 
respectively.  The  specific  case  of  p  =  2  is  the  astroid  curve  shown  in  (a)  and  investigated 
previously  for  backscattering  from  freely  rising  bubbles  in  water.4  (See  also  Sec.  HI  of  the 
present  report)  The  orientation  of  the  x  axis  is  parallel  to  the  direction  \|/  =  0  in  the  polar 
coordinate  system  used  to  specify  W.  The  number  of  (nondegenerate)  cusp  points  for  odd 
values  of  p  is  fewer  than  might  have  been  anticipated.  For  example,  while  p  =  2 
produces  4  cusp  points,  Fig.  3(a),  the  caustic  for  p  =  3,  Fig.  3(d)  has  only  3  cusp  points. 

Amott12*13  also  calculated  the  acoustic  wavefield  which  decorates  the  caustic  for  the 
specific  case  of  p  =  2.  To  display  the  wavefield  the  squared  modulus  of  the  pressure  is 
calculated  and  imaged  with  a  grey  scale  where  bright  corresponds  to  large  acoustic 
intensity.  Figure  4  shows  the  resulting  acoustic  wavefield  patterns  for  parameters  a/b  = 
0.25,  5/b  =  -0.02,  and  (a)  kb  =  100,  (b)  kb  =  500,  and  (c)  kb  =  1000  where  k  =  2%fk. 
Superposed  on  these  patterns  is  a  white  curve  which  marks  the  astroid  caustic.  For 
comparison,  an  acoustic  unperturbed  toroidal  wavefront  produces  a  pattern  of  concentric 
circular  rings  since  the  pressure  amplitude  is  proportional  to  J0(kb  siny)  where  y  is  the 
angle  tan-1[(x2  +  y2)l/2/z],  (See,  e.g.,  the  glory  scattering  pattern  of  Lamb  waves  on  a 
spherical  shell.4*5’8-9)  The  wavefields  are  changed  as  a  consequence  of  the  perturbation 
f(\jr)  of  the  outgoing  wavefront  The  pattern  in  Fig.  4(c)  near  the  cusp  point  has  features 
like  that  of  a  Pearcey  pattern  discussed  previously  and  later  in  Sec.  n. 

C.  Observation  of  Lips  Caustics  in  Light  Backscattered  from  Oblate  Drops 

Catastrophe  theory  allows  the  caustics  leaving  two  cusp  points  to  join  smoothly.15 
The  resulting  bounded  caustics  have  the  appearance  of  a  pair  of  lips.  (For  comparison, 

Fig.  3(a)  and  (d)  show  the  smooth  joining  of  caustics  from  4  and  3  cusp  points, 
respectively.)  With  an  additional  flattening  of  the  outgoing  wavefront,  the  lips-like  caustics 
are  drawn  together  and  the  two  cusp  points  gradually  merge.  The  merging  of  cusp  points 
is  known  as  a  "lips  event"  in  catastrophe  terminology.15 
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Nye15  predicted  that  lips  caustics  occur  in  light  backscattered  from  oblate  liquid 
drops.  Since  manifestations  of  lips  caustics  are  relatively  unexplored  and  the 
backscattering  of  light  or  sound  from  penetrable  spheroids  is  a  problem  of  general  interest, 
an  experimental  investigation  of  these  phenomena  was  undertaken.  The  initial  research  was 
the  M.S.  degree  project  of  Harry  Simpson16  though  it  is  anticipated  that  the  scope  will  be 
expanded  to  give  a  Ph.D.  project.  All  experiments  to  date  have  been  carried  out  by 
observing  light  backscattered  from  oblate  acoustically  levitated  drops  of  water  for  which  the 
refractive  index  |i  =  1.332.  The  first  task  was  to  build  a  new  acoustic  levitator  which 
would  allow  a  beam  splitter  to  be  placed  close  to  the  vibrating  surfaces  of  the  levitator.  The 
apparatus  is  shown  in  Fig.  5.  The  summary  of  the  work  below  is  taken  from  a  recent 
overview  of  research  on  optical  catastrophes.17 

Nye  predicted  two  lips  events  occur  in  light  backscattered  from  spheroidal  drops 
having  critical  aspect  ratio  D/H  of5 

qL1  *  [m/C 2p-2)]1/2-  1.416,  qL2  =  [(2n-l)/(2^t-2)] ^  »  1.584  (18,19) 

where  D  is  the  diameter  of  the  drop  and  H  is  it  height.  (See,  e.g..  Fig.  3a  of  Ref.  18.) 

For  qLi  <  q  <  qL2.  where  q  =  D/H,  there  are  no  caustics  in  the  near  backward  direction. 
Marston's  analysis  of  the  cusp  point  location18  predicts  that  the  horizontal  scattering  angle 
03  to  the  cusp  point -4  180°  as  q-»qu-  That  analysis  indicates  e.g.  that  if  q  is 
reduced  below  qLi  to  1.40,  then  103  -  180°l  =*  26°. 

Nye's  derivations  are  brief  so  a  discussion  is  merited.  The  condition  for  the  event 
at  qLj  is  that  paraxial  rays  in  the  vertical  plane  of  symmetry  are  focused  on  the  back  of  the 
drop  (Fig.  6).  The  focal  length  in  this  plane  is  given  by  the  lens  maker's  formula  f  = 
Mp/(p.-l)  where  p  =  H2/2D  is  the  vertical  radius  of  curvature  at  the  equator.  Setting  f  = 

D  and  solving  for  D/H  yields  Eq.  (18).  The  condition  for  the  event  at  <\\ji  is  shown  in 
Fig.  7.  The  internal  paraxial  focus  F  is  now  a  distance  p  from  the  back  side  of  the  drop 
so  that  the  reflected  rays  are  also  focused  at  F.  Setting  p  =  D-f  yields  Eq.  (19). 


Fig.  5.  Top  view  of  apparatus.  Light  is  backscattered  from  the  drop  through  a  beam 
splitter  into  a  camera  which  records  the  backscattering  pattern. 


Fig.  6.  Profile  of  oblate  drop  in  a  vertical 
plane  containing  the  vertical  axis  of  rotational 
symmetry.  Rays  are  focused  at  the  back  of 
the  drop  for  the  1st  lips  event 


Fig.  7.  Profile  as  in  Fig.  6  but 
showing  the  focusing  of  rays 
for  the  2nd  lips  event 
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We  have  confirmed  these  events  exist  by  viewing  laser  light  backscattered  from 
levitated  drops  of  water.  ^  The  ultrasonic  field  for  our  levitator  is  excited  by  a  pair  of 
vibrating  rods  which  face  each  other.  This  facilitates  mounting  a  camera  (focused  on 
infinity)  closer  to  the  drop  than  previously  possible.  The  incident  beam  reflects  from  a 
beam  splitter  (mounted  adjacent  to  the  rods)  onto  the  drop.  The  backscattered  light  passes 
through  the  splitter  to  the  camera.  A  second  camera  records  the  drop's  profile.  Figure  8 
shows  backscattering  patterns  recorded  for  various  D/H.  (The  bright  point  visible  in  these 
is  light  from  a  distant  beam  dump  illuminated  by  that  part  of  the  beam  not  scattered  by  the 
drop.  It  is  useful  as  a  fiducial  mark  of  the  backward  direction.)  Drops  flattened  by  this 
amount  from  acoustic  radiation  pressure  may  have  significant  deviations  from  the  perfect 
spheroidal  shape  assumed  in  the  derivations.  Hence  it  is  appropriate  to  consider  an 
optically  equivalent  spheroid  (with  aspect  ratio  q^  where  the  vertical  curvature  at  the 
equator  p-1  matches  that  of  the  actual  drop.  Let  qsj  denote  the  apparent  qs  for  the  jth 
part  of  Fig.  8.  The  patterns  in  Fig.  8  may  be  understood  as  follows:  qsa  is  sufficiently 
less  than  qu  that  the  lips  appear  as  two  horizontal  caustic  lines.  These  are  the  distinct 
bright  lines  which  are  displaced  in  (a)  by  vertical  angle  -1.2°.  The  drop  is  flatter  in  (b)  so 
that  qS5  is  much  closer  to  qu  and  the  horizontal  caustics  are  nondistinct.  Nevertheless, 
as  in  (a),  03  is  not  close  enough  to  180°  for  the  cusp  points  (where  the  caustics  terminate) 
to  be  in  the  field  of  view.  (A  somewhat  wider  field  of  view  than  the  present  17°  would  be 
needed  to  photograph  symmetric  cusp  points.)  In  (c),  a  slight  increase  of  the  acoustic 
amplitude  raises  qs  to  q#  >  qu  so  that  no  caustics  are  visible.  The  reason  why  the  actual 
merging  of  symmetric  cusp  points  at  180°  would  be  difficult  to  photograph  is  evident  by 
inspection  of  Fig.  4  of  Ref.  18:  the  rate  at  which  63  -*  180°  diverges  as  qs  -»  qu-  An 
additional  flattening  of  the  drop  gives  rise  to  a  new  horizontal  caustic  at  qsd  >  qL2  because 
of  the  second  lips  event. 

To  facilitate  a  quantitative  check  of  the  theory,  several  drops  were  photographed  for 
patterns  equivalent  to  (a)  -  (d)  of  Fig.  8.  For  each  photo,  qs  was  determined  by  measuring 
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the  equatorial  radius  of  curvature  p  and  averaged  with  other  measurements  of  the  jth 
class.  The  resulting  averages  are  qsa  =  1-398  ±  0.085,  qsb  =  1.408  ±  0.052,  q^  =  1.461 

±  0.048,  and  q^  =  1.587  ±  0.093.  These  results  are  consistent  with  the  predictions  that 
Qsa  <  qu>  qsb  5  qu>  qu  5  qsc  <  qL2  and  qsd  >  qLi.  The  size  of  the  uncertainties  is 

indicative  of  the  difficulty  in  determining  p  from  photographs.  The  D  ranged  from  1.3 
to  2.3  mm.  In  spite  of  the  aforementioned  uncertainties,  the  predicted  sequence  of  lips 
caustics  is  confirmed.  This  is  the  first  observation  of  these  events  for  scattering  from 
oblate  objects. 

D.  Observations  of  Acoustical  and  Optical  Transverse  Cusps  Produced  by 
Reflection 

In  the  previous  Annual  Report,4  experiments  are  described  in  which  a  burst  of 
ultrasound,  only  a  few  cycles  in  duration,  was  reflected  from  a  curved  polished  metal 
surface  in  water.  The  surface  was  curved  in  the  generic  shape  to  produce  transverse  cusp 
caustic.14-19  See.  Fig.  9.  Depending  on  the  location  of  the  receipting  hydrophone,  distinct 
echoes  could  be  resolved  in  time;  echoes  merged  when  the  caustic  was  crossed. 

The  research  during  the  contract  period  under  consideration  was  to  explore  the 
reflected  wavefield  for  a  quasi-steady-state  incident  burst  of  ultrasound.  This  work  was 
part  of  the  M.S.  degree  project  of  Carl  Frederickson,20  though  it  is  anticipated  that  the 
scope  will  be  expanded  to  that  of  a  Ph.D.  project  The  use  of  long  tone  bursts  has  two 
advantages  over  the  work  reported  in  Ref.  4.  (i)  A  theory  for  the  steady  state  reflected 
wavefield  was  given  by  Marston4-19  which  predicts  the  pattern  is  that  of  a  Pearcey 
function14 


P(X,Y)  =  j*  exp[i(s4  +  s2X  +  sY)]ds, 


(20) 
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Fig.  9.  The  wave  from  a  point  source  reflects  to  give  a  transverse  cusp  in  the  observation 
plane.  The  reflecting  surface  illustrated  has  a  height  relative  to  the  xy  plane  of  h(x,y)  = 
cix2  +  C2y2x  with  ci  <  0  and  C2  <  0. 


Fig.  10.  (on  the  next  page.)  Observed  and  predicted  transverse  cusp  caustics  and 
wavefield  patterns  for  a  curved  polished  metal  surface  in  water  (all  on  the  same  scales),  (a) 
Shows  the  optical  caustic  recorded  by  placing  a  light  source  in  the  water  and  scanning  a 
photocell,  (b)  Shows  the  acoustical  transverse  cusp  pattern  for  the  same  reflector  where 
sound  frequency  is  610  kHz.  (c)  Shows  the  calculated  pattern  for  the  acoustic  wavefield 
for  that  frequency  of  sound  (see  text).  The  horizontal  width  of  the  region  displayed  in  (a) 
is  roughly  9  cm  while  the  vertical  width  is  roughly  14  cm;  zs  =  133  cm  while  z  *>  48  cm. 
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where  X  and  Y  are  related  to  the  transverse  coordinates  (u,v)  of  the  observer  (Fig.  9). 

(ii)  For  the  laboratory  scale  experiments  described  in  Ref.  4,  it  was  difficult  to  fully  resolve 
distinct  echoes  because  of  the  bandwidth  restrictions  of  the  transducer  and  amplifier 
systems  used.  This  resolution  problem  is  avoided  when  studying  the  structure  of  the 
acoustic  caustic  by  the  present  method  since  the  incident  burst  is  selected  to  be  long  enough 
(burst  duration  %  >  400  |xs)  that  all  of  the  echoes  overlap  completely  in  the  region  of  time 
and  space  sampled  for  the  reflected  wavefield.  (The  present  method  is  analogous  to 
infening  small  path  length  differences  with  an  interferometer  instead  of  by  direct 
measurement  in  the  time  domain.)  Marston's  theory4*19  relates  X  and  Y  in  Eq.  (20)  not 
only  to  (u,v,z)  of  the  observer  but  also  to  the  source  coordinates  (us,  vs,  zs),  the 
wavelength  X,  and  the  parameters  Cj  of  the  elevation  function  of  the  reflecting  surface 

h(x,y)  =  CjX2  +  c2y2x  +  c3y2.  (21) 

The  wavefronts  leaving  the  source  at  (us,  vs,  zs)  are  approximately  spherical  for  those 
waves  diverted  toward  the  surface. 

In  the  experiment,  the  receiver  hydrophone  was  scanned  in  a  raster  pattern  with  the 
aid  of  a  modified  "x,y"  plotter.  The  received  signal  was  amplified  and  rectified  and 
sampled  at  an  appropriate  time.  The  resulting  sample-and-hold  voltage  produced  was 
proportional  to  the  pressure  amplitude  at  the  current  hydrophone  position.  This  voltage 
was  used  to  modulate  the  intensity  of  an  oscilloscope  beam  subjected  to  a  raster  scan  in 
synchronization  with  the  hydrophone.  Figure  10(b)  shows  a  representative  photograph  of 
the  pattern  on  the  oscilloscope  which  images  the  sound  field.  Bright  regions  correspond  to 
the  regions  of  greatest  acoustic  pressure. 

The  reflecting  surface  used  was  fabricated  with  ci  *  0  and  C3  *  0  in  Eq.  (21). 
One  approach  to  testing  the  theory  would  require  accurate  measurement  of  C2,  us,  and  vs. 
We  chose  a  different  method  which  involves  comparison  with  the  optical  caustic  produced 
by  placing  an  optical  point  source  at  the  same  location  as  the  acoustic  source.  This  was 
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done  with  the  aid  of  an  optical  fiber  (driven  by  a  laser)  with  its  output  end  in  the  water  tank. 
The  receiver  hydrophone  was  replaced  by  a  water-proof  photo  detector.  The  output  of  the 
detector  modulates  the  oscilloscope.  This  gives  the  image  of  the  optical  caustic  shown  in 
Fig.  10(a).  It  should  be  emphasized  that  all  of  the  geometrical  properties  are  the  same  in 
Fig.  10(a)  and  (b);  the  difference  in  the  patterns  is  primarily  a  consequence  of  X  differing 
by  a  factor  of 5200.  Let  Mu  and  Mv  denote  the  horizontal  and  vertical  magnification 
factors  for  these  patterns  (relative  to  actual  laboratory  coordinates).  For  these  patterns,  Mu 
<  Mv  while  the  respective  values  of  Mu  and  Mv  are  the  same  in  Fig.  10(a),  (b),  and  (c). 
The  cusp  curve  of  the  caustic  is  characterized  by4*19  an  "opening  rate."  By  overlaying  such 
curves  on  Fig.  10(a),  an  empirical  opening  rate  was  determined  where  use  is  made  of  Mu 
and  Mv.  Then  Mansion's  theory  was  used  to  give  the  surface  parameter  C2  in  Eq.  (21) 
and  the  acoustic  wavefield  shown  in  Fig.  10(c).  The  agreement  between  Figs.  10(b)  and 
(c)  is  satisfactory,  though  experiments  are  under  way  to  repeat  this  comparison  over  a 
range  of  frequencies.  For  smaller  acoustic  X,  the  structure  in  the  pattern  is  more  closely 
spaced.  For  one  of  the  reflecting  surfaces  C2  could  be  accurately  estimated  directly  and 
this  gave  a  C2  in  satisfactory  agreement  with  the  one  determined  optically. 

To  evaluate  the  theory,  it  was  necessary  to  use  a  numerical  algorithm  to  evaluate  the 
Pearcey  function,  Eq.  (20).  Another  graduate  student  supported  by  this  contract,  Cleon 
dean,  participated  in  that  accomplishment. 

III.  OPTICS  OF  BUBBLES  IN  WATER 

A.  Physical  Optics  of  Bubbles  in  Water 

Previous  Annual  Reports4-21  have  surveyed  research  into  the  physical  optics  of 
bubbles  in  water.  Contrary  to  what  may  have  been  anticipated  there  can  be  significant 
errors  if  purely  geometrical  methods  (adapted  to  retain  phase  information)  are  used  to 
calculate  the  scattered  irradiances.  This  is  because  even  though  bubbles  are  very  much 
larger  than  the  wavelength  of  light  in  water  Xw,  the  effects  of  diffraction  are  significant 
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near  caustics  and  near  the  critical  angle.  Furthermore,  surface  coatings  which  can  be 
present  on  bubbles  in  ocean  water,  may  affect  the  scattering  in  ways  which  merit 
exploration.4*21*22 

The  emphasis  on  this  section  will  be  on  accomplishments  during  the  past  year 
which  are  yet  to  be  published.  Marston  was  recently  asked  to  give  an  overview  of  bubble 
optics23  and  parts  of  the  summary  below  are  adapted  from  that  manuscript. 

B.  Brewster  Angle  Scattering  from  Bubbles:  Theory  and  Experiment 

When  the  electric  field  is  polarized  parallel  to  the  plane  of  incidence,  the  Fresnel 
reflection  coefficient  of  a  flat  clean  interface  vanishes  when  the  angle  of  incidence  is  at 
Brewster's  condition22 

iB  =  arctan(n./nw),  Brewster  scattering  angle  =  0B  =  180°  -  2ig  =  106.2°  (22,23) 

where  ni  and  nw  are  the  refractive  indices  of  air  and  water.  Equation  (23)  follows  from 
geometry  and  the  assumption  that  Brewster's  condition  is  not  significantly  affected  by 
curvature  provided  the  bubble  radius  a  »  X.w.  las  local  reduction  in  the  reflectivity 
manifests  itself  in  calculated  I2,  the  normalized  irradiance  for  this  polarization  for  an 
uncoated  bubble.  This  was  shown  in  Fig.  9  of  last  year’s  Annual  Report.4  The  visibility 
of  fine  structure  fringes  is  greatly  reduced  when  the  scattering  angle  0  is  close  to  106.2°. 
Similar  calculations  done  for  coated  bubbles  predict,  however,  a  much  smaller  reduction  in 
the  fringe  visibility  for  0  near  106.2°.  Consequently,  if  Brewster  angle  effects  can  be 
observed  for  freely  rising  uncoated  bubbles,  it  is  plausible  that  such  effects  could  be  used 
to  optically  discriminate  uncoated  from  coated  bubbles.  Stefan  Baumer  carried  out  the  first 
observations  of  Brewster  angle  light  scattering  from  air  bubbles  rising  in  water  with  the 
support  of  this  contract.  His  M.S.  thesis  which  describes  these  experiments  has  been 
issued  as  a  Technical  Report.24 

The  scattering  patterns  in  the  angular  region  near  106°  were  photographed  for 
bubbles  rising  freely  through  distilled  water.  This  was  done  to  verify  for  the  first  time,  that 
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a  local  reduction  in  fringe  visibility  is  observable.  The  bubbles  were  illuminated  by  a 
horizontally  propagating  300  mW  beam  from  an  Ar-ion  laser  for  which  the  wavelength  in 
air  Xj  =  514.5  nm  and  k  =  2jtnw/(ni  514  nm) «  16.29  pur1.  The  already  polarized  output 
of  the  laser  was  made  more  perfectly  polarized  (with  a  horizontal  E  field)  by  passing  the 
beam  through  an  ellipsometer-grade  polarizer.  Light  in  the  horizontal  scattering  plane 
passed  through  a  window  into  a  Polaroid  filter  and  camera  focused  on  infinity.  The 
apparatus  was  generally  similar  to  the  ones  used  to  view  critical  angle  scattering  except  for 
changes  needed  to  facilitate  viewing  the  Brewster  region  and  the  reduction  in  the  level  of 
spurious  and  background  scattering.  After  rising  through  the  beam,  the  bubble  floated 
against  a  glass  slide  and  its  radius  was  measured  directly.  The  photographic  negatives  of 
the  scattering  were  scanned  with  a  microdensitometer  from  which  an  exposure  En  as  a 
function  of  the  scattering  angle  6  could  be  determined.  The  angle  calibration  of  the 
photograph  was  achieved  by  use  of  a  mirror-goniometer  technique.  Even  with  no  bubble 
in  the  beam,  the  photograph  would  be  weakly  exposed.  The  smoothed  irradiance  pattern 
inferred  from  a  background  photograph  will  be  denoted  by  f(0).  It  was  desired  to 
compare  the  scattering  patterns  from  bubbles  with  predictions  from  Mie  theory  for  clean 
(uncoated)  bubbles.  To  facilitate  this  comparison,  it  was  necessary  to  introduce  three 
positive  adjustable  parameters  a,p,  and  C  from  which  the  normalized  irradiance  12(0) 
could  be  inferred  from  the  exposure  En(0)  for  a  given  bubble  via 

I2(0)  =  P[en(0)  -  af(0)  +  C].  (24) 

The  parameter  a  was  selected  so  as  to  remove  a  component  of  En  attributable  to 
background;  C  was  selected  to  eliminate  an  apparent  offset  The  scaling  constant  P  was 
adjusted  to  facilitate  comparison  with  a  locally  averaged  form  of  12(0)  from  Mie  theory.  A 
similar  local  average  was  performed  on  the  "raw"  En(0)  from  the  negatives  so  as  to 
reduce  noise  associated  with  the  graininess  of  each  photograph. 
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Fig.  11.  Normalized  irradiance  for  a  bubble  of  radius  78.5  ^im  in  green  light:  (a)  is 
inferred  from  a  photograph  and  confirms  the  local  reduction  of  fringe  visibility  near  the 
Brewster  condition;  (b)  is  a  narrow  angular  average  if  I2  from  Mie  theory  for  ka  =  1278. 
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Figure  1 1  is  a  representative  comparison  of  data  with  theory  for  a  case  where  the  local 
average  was  over  ±  0.04°.  Outside  of  the  region  from  about  105°  to  108°,  the  quasiperiod  of  the 
irradiance  oscillations  is  fairly  uniform.  Measurement  of  that  quasi-period  shows  that  the  oscil¬ 
lations  are  primarily  associated  with  the  interference  of  the  surface  reflected  ray  with  a  (far-side) 
2  chord  ray.  For  0  close  to  0b  =  106°,  there  is  a  noticeable  reduction  in  the  fringe  visibility 
and  the  quasi-period  changes.  The  reflected  ray  is  suppressed  there  and  the  fringes  can  be 
attributed  to  the  interference  of  other  (weak)  rays.24  The  important  result  is  that  the  same  general 
features  are  observable.  It  is  noteworthy  that  the  local  reduction  of  the  fringe  visibility  near  0  b 
is  evident  in  the  "raw"  En(0).  Similar  detailed  comparisons  were  earned  out  for  bubbles  having 
radii  ranging  from  60  to  80  4m,  while  reduction  in  fringe  visibility  was  noted  with  radii  from  50 
to  100  pm.  Since  the  experiments  were  carried  out  in  distilled  water,  no  optically  significant 
coating  was  expected  and  no  evidence  for  such  a  coating  was  found.  This  confirmation  of 
Brewster  effects  for  uncoated  bubbles  suggest  that  the  effects  of  a  film  may  be  observable. 

Since  geometric  features  are  important  near  the  critical  scattering  angle  0C  =  82.8° 
and  0  b,  it  may  be  argued  that  the  ratio  <I2(9c)>/<I2(8b)>  depends  only  weakly  on  ka 
for  bubbles  where  <  >  denotes  a  local  average  about  the  indicated  scattering  angle. 
Measurement  of  that  ratio  may  be  useful  for  discrimination  of  bubbles  from  particles-in- 
water.26 

C.  Backscattering  from  Freely  Rising  Spherical  and  Spheroidal  Air 
Bubbles  in  Water:  Glory  Scattering  and  the  Astroid  Caustic 

For  a  spherical  bubble,  some  of  the  rays  incident  with  non-vanishing  impact 
parameters  leave  the  bubble  scattered  exactly  backwards.26  As  a  consequence  of  the 
azimuthal  symmetry  of  spheres,  the  scattering  associated  with  such  rays  is  weakly  focused 
along  the  backward  axis.  This  weak  focusing  is  generally  termed  "axial  focusing"  or 
"glory  scattering"  and  manifests  itself  in  a  local  enhancement  of  backscattering  from 
spherical  cloud  droplets;  the  detailed  cause  of  glory  rays  is  different  for  drops  and  bubbles. 
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The  shape  of  the  backward  directed  wavefront  is  toroidal  and  is  of  the  form  described  by 
W(s)  in  Eq.  (13). 

Real  bubbles  rising  through  water  take  on  an  oblate  shape  as  a  consequence  of 
hydrodynamic  forces  on  the  bubble.  For  a  sufficiently  oblate  bubble,  the  outgoing 
wavefront  becomes  distorted  as  shown  in  Fig.  9  of  a  previous  Annual  Report21  and  the 
focal  properties  of  the  scattering  are  changed.  The  wavefront  shape  becomes  perturbed 
such  that  it  is  locally  well  approximated12*13  by  Eqs.  (13)  -  (15)  with  p  =  2  in  Eq.  (14). 
Hence  the  axial  caustic  unfolds  to  give  an  astroid  caustic  having  the  shape  shown  in  Fig. 
3(a).  Because  of  polarization  effects  and  the  superposition  of  fields  from  several  such 
wavefronts,  the  scattering  patterns  are  much  more  complicated  than  those  shown  in  Fig.  4 
for  the  acoustic  case. 

During  the  contract  period  under  consideration  Amott  photographed  and  modeled 
backscattering  patterns  from  bubbles  rising  through  water  which  were  sufficiently  oblate  to 
display  an  astroid  caustic.  The  research  is  described  in  his  Ph.D.  dissertation12  and  in  a 
short  publication27  and  will  be  only  briefly  summarized  here.  All  of  his  observations  are  of 
cross-polarized  backscattering  patterns.  In  the  ones  described  in  previous  Annual 
Reports4*21  the  incident  light  was  vertically  polarized  so  that  the  polarization  is  symmetric 
relative  to  the  rotational  axis  of  the  oblate  bubble  (since  the  latter  axis  is  also  vertical).  In 
the  new  observations,  the  incident  light  was  polarized  at  45°  relative  to  the  vertical  which  is 
asymmetric  relative  to  the  bubble.  With  this  change  of  illumination,  the  backscattering 
pattern  is  more  intense  near  the  cusp  points  of  the  astroid  than  in  the  previous  observations. 
Amott  extended  his  previous  calculation  of  the  relevant  diffraction  integrals  to  include  the 
case  of  asymmetric  polarization.  Figure  12  shows  the  observed  (a)  and  calculated  (b) 
pattern  for  asymmetric  illumination  of  a  bubble  having  a  diameter  d  =  654  pm.  The 
caustics  associated  with  the  3  and  4  chord  glory  wavefronts  are  highlighted  in  (b).  The 
theory  reproduces  the  general  features  and  dimensions  of  the  observed  pattern  if  not  its 
detail. 


Rg- 12.  Observed  (a)  and  calculated  (b)  backscattering  patterns  for  an  oblate  bubble  rising 
vertically  through  water  with  a  diameter  d  -  654  m.  The  vettical  rise  dintcrion  is  shown 

by  the  arrow.  The  illumination  was  green  laser  light  polariaed  at  45  deg  teladve  to  the 
vertical.  Astroid  caustics  are  superposed  on  (b). 
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In  related  work,  Amott  measured  the  apparent  caustic  dimensions  for  patterns  with 
vertically  polarized  illuminations.  He  used  those  measurements  to  infer  the  oblateness  of 
the  bubble  where  use  is  made  of  the  optical  model.  The  inferred  oblateness  is  in  fair 
agreement  with  the  results  of  a  hydrodynamics  approximation.  I2*27  This  method  of 
determining  the  oblateness  of  a  rising  bubble  is  very  much  more  sensitive  than  simple 
photography  of  a  rising  bubble.  The  difference  between  the  major  and  minor  diameters 
may  be  inferred  even  when  that  difference  is  <  1%. 

D.  Asymptotic  Series  for  Scattering  at  the  Critical  Angle  from  CAM 
Theory 

While  the  physical  optics  approximation28  (POA)  gives  a  simple  and  quantitative 
understanding  of  coarse  structure  present  is  scattering  with  0  <  8C  =  82.8°,  the  critical 
value  for  total  refection,  comparisons  with  Mie  theory29  show  that  the  POA  noticeably 
underestimates  the  inradiance  at  0C  when  ka  is  quite  large  (roughly  ka  >  2000).  The 
error  is  a  consequence  of  the  auxiliary  assumptions  made  in  the  POA  so  that  the  relevant 
diffraction  integral  could  be  reduced  to  the  form  of  a  Fresnel  integral.28  To  facilitate  a  more 
rigorous  analysis  Ferrari  and  Nussenzveig30*31  extended  the  modified  Watson 
transformation  of  the  Mie  series  to  the  case  of  bubbles.  The  results  is  known  as  the 
complex  angular  momentum  (CAM)  theory.  They  obtained  excellent  agreement  with 
angular  scattering  patterns  computed  from  Mie  theory.30  The  important  result  is  an 
approximation  for  the  amplitude  of  reflected  light  in  terms  of  integrals  denoted  as  Pp  and 
Fp  and  given  the  names  Pearcey-Fock  and  Fresnel-Fock  integrals.  The  integrals  depend 
on  ka,  0,  and  0C.  The  integrals  were  evaluated  numerically  by  Ferrari.30 

To  obtain  a  simple  series  approximation  to  the  scattering  at  0  =  0C,  Cleon  Dean 
found  series  expansions  of  Pp  and  Ff  in  the  special  case  0  =  0c-  This  work  was  only 
partially  supported  by  this  contract32  The  contribution  of  the  reflected  light  to  the 
normalized  inradiance  becomes 
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q=0  q=0 

where  0  =  ka,  Aj,o  =  Bjfo  =  1/2  and  the  coefficients  Aj,q  and  Bj,q  are  complex  for  q  >  1 
and  j  is  a  polarization  index.  Note  that  0*  ^  =  1  and  0*  ^  =  1  for  q  =  0.  Because  the 
analysis  leading  to  the  Aj>q  and  Bj,q  is  tedious,  it  was  necessary  to  terminate  the  series  at 
small  Ma  and  Mb-  While  the  form  of  the  series  is  certain,  some  aspects  of  the  analysis 
for  the  Aj>q  and  Bjfq  with  q  >  0  are  being  reexamined  as  of  the  time  of  this  writing. 

Some  of  the  noteworthy  features  evident  from  inspection  of  Eq.  (25)  are  (i)  lj(0c,  0  — > 

— >  1  which  is  the  correct  geometric  optics  limit  for  total  reflection  and  (ii)  the  leading 
correction  to  the  geometric  optics  limit  is  0(0* 1/6).  This  is  of  fundamental  interest  since 
away  from  special  angles,  the  corrections  to  geometrical  optics  are  usually  0(0*1).  Hence 
0  must  be  usually  large  for  0  =  0C  for  Ij  to  be  close  to  1.  For  example,  for  j  =  1, 
which  is  E  field  parallel  to  the  scattering  plane,  diffraction  is  especially  important  and 
<l2>  from  Mie  theory  has  only  reached  0.62  for  ka  =  20000,  corresponding  to  a  bubble 
radius  a  of  1.5  mm  in  red  light. 

This  analysis  constitutes  a  part  of  Cleon  Dean's  dissertation  project  which  should 
be  completed  in  the  next  few  months.  In  addition  to  being  a  fundamental  problem  in 
scattering  physics  because  0(0*  ^)  terms  vanish  slowly,  the  analysis  has  applications.23 
It  is  important  for  a  full  understanding  of  the  unique  monotone  increasing  relationship 
between  scattered  power  and  bubble  radius  for  0  near  0C.  That  relation  has  applications 
in  instruments  to  measure  bubble  size  and  pulsation  amplitudes.  Furthermore,  there  is  a 
close  analogy  with  critical  angle  reflection  of  sound  from  curved  fluid  interfaces,33  so  the 
analysis  could  be  extended  to  give  an  improved  understanding  of  acoustic  reflections  if 
such  is  merited. 
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With  the  support  of  this  contract  C.  Dean  has  also  been  working  on  applying 
wavefront  tracing  to  a  situation  in  catastrophe  optics  like  that  described  in  Ref.  1 8.  There  is 
no  significant  new  progress  on  that  problem. 

IV.  ACOUSTICAL  PHASE  CONJUGATION  AND  ASPECTS  OF 
ACOUSTIC  WAVEFRONT  REVERSAL 

A.  Theory  for  Wavefront  Reversal  Via  Three-Wave  Mixing  in  a  Layer  of 
Enhanced  Nonlinearity 

Three-wave  mixing  is  usually  not  the  method  of  choice  in  optical  phase  conjugation 
as  a  consequence  of  dispersion.34.  It  is,  however,  potentially  useful  as  a  method  of 
acoustical  phase  conjugation  where  a  layer  of  bubbles  is  used  to  enhance  the  region  of 
nonlinearity.35-36  To  improve  our  understanding  of  focal  point  locations  and  the  relevant 
experimental  parameters,  an  analysis  of  the  difference  frequency  signal  produced  is  carried 
out.  This  is  summarized  below.  A  brief  overview  appears  at  the  end  of  this  section. 

The  geometry  considered  is  shown  in  Fig.  13.  Let  Si  denote  the  location  of  the 
pump  wave  of  radian  frequency  coi,  and  let  S2  denote  the  location  of  the  probe  wave  of 
radian  frequency  ©2.  It  is  desired  to  produce  a  wave  which  converges  back  generally 
towards  S2  by  way  of  three-wave  mixing.  The  frequency  of  the  reversed  wave  is  ©3  = 
©l  -  002-  The  Cartesian  coordinates  of  the  jth  source  will  be  denoted  by  (xj,  yj,  zj),  j  =  i 
and  2  and  we  may  take  S2  at  (0, 0,  Z2)  and  Si  at(xi,0,  zi).  The  pressure  due  to  these 
sources  is: 

p(x,y,z,t)  =  ■i[Ple<°1  +P2e°>z  +  c.c.],  (26) 

where  c.c.  denotes  the  complex  conjugate  and  Pj  =  Aj  exp(-ikjRj)  where  kj  =  (Oj/c,  Aj  is 
the  source  strength,  and  Rj  is  the  distance  of  (x,y,z)  from  the  Si  or  S2. 
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Layer  of  Microbubbles 


Fig.  14.  Current  experimental  configuration 
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The  present  calculation  considers  parametric  sources  of  sound  where  the  source 
strength  density  q  is  associated  with  a  quadratic  nonlinearity.  Westervelt's  parametric 
source  equation  reduces  to37 

vV-ip^-PoV  (27> 

c 


q  =  (p  c4)'1  e  3(p2)/3t 


(28) 


where  the  subscript  t  denotes  differentiation  with  respect  to  time  and  hydrodynamic 
corrections  are  neglected  so  that  Eq.  (28)  describes  the  distributed  source  associated  with 
the  quadratic  nonlinearity  in  the  equation  of  state.  (In  conventional  notation,  the 
nonlinearity  parameter  e  becomes  B/2A  for  a  homogeneous  fluid.) 

Inspection  of  Eqs.  (26)-(28)  shows  that  a  component  of  parametrically  generated 
wave  p'  oscillates  at  003  with  q  being  proportional  to  P2*(x,y,z)  Pi(x,y,z).  The 
relevant  component  of  p'(x',y',z\t)  is  given  by  the  real  part  of 
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<$(r,t)  =  (c^)  P^  e  e(©  ,co2,r)/poc  .  (30) 

The  integral  in  (29)  extends  over  the  region  of  the  parametric  source.  For  a  bubbly  liquid, 
e  becomes  complex  to  take  into  account  the  phase  of  the  response  of  the  bubbles.37’38 
Furthermore  lei  may  be  »  B/2A  of  the  surrounding  water  and  for  the  purposes  of 
estimating  the  amplitude  of  the  reversed  wave  produced  in  a  bubble  layer  we  may  take  e  = 

0  outside  of  the  layer.  We  will  be  interested  in  P3  generated  by  the  response  of  bubbles 
in  a  flat  thin  layer  of  thickness  L.  Let  the  layer  lie  in  the  plane  z  =  zo.  For  sufficiently  thin 
L, 

efw,  to2,  r)  =  e(co  ,  0^)  L  8(z  -  zo), 


(31) 
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where  8(z  -  Zq)  is  a  Dirac  8  function  so  that  (29)  reduces  to  a  two-dimensional  integral. 
That  P3  describes  a  reversed  wave  may  be  seen  by  use  of  the  parabolic  (or  Fresnel) 
approximation  to  simplify  the  various  propagation  distances  which  arise  (such  as  the 
distances  Rj  from  the  jth  source  to  a  given  region  of  the  bubble  layer.  For  simplicity  we 
take  the  plane  of  the  layer  to  be  at  Zq  =  0  and  Eqs.  (29)  -  (31)  reduce  to 

-(o^Le(co  co  )A, 

P3(x\y',z’)  =  -2 - L-l_i  J(x',y',z')  (32) 
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where  X2  =  0  for  the  case  considered.  Inspection  of  (33)  shows  that  J  diverges  if  F  =  0, 
G  =  0,  and  H  =  0.  These  conditions  may  be  used  to  directly  specify  the  (x',y',z')  where 
the  reversed  wave  comes  to  a  focus  according  to  geometrical  optics.  The  focus  is  at  y'  =  0 
with 
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X  =  -  X, 
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22  1  r 


Notice  that  when  z\  -»  °°  and  xj  =  0  (so  that  the  pump  wave  is  a  plane  wave  at  normal 
incidence)  and  coi  =  2c02  (so  that  C03  =  C02)  then  z'  =  Z2  and  x'  =  0.  The  reversed  wave 


is  focused  back  on  S2  which  is  a  requirement  of  true  phase  conjugation.  Specific  cases  of 
Eq.  (35)  for  the  "focal  point  shifts"  had  been  derived  previously  by  Marston4  by  less 
transparent  methods.  Equation  (35)  also  agrees  with  results  in  Ref.  36. 

A  mathematical  demonstration  of  the  conjugate  nature  of  the  reversed  wave  in  the 
aforementioned  special  case  of  ©2  =  coi/2  =  kc,  z\  -» «>  with  xj  =  0,  follows  by  the 


stationary  phase  approximation  (SPA)  of  J  in  Eq.  (32)  which  gives 
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where  R  is  the  distance  of  the  observation  point  (x',y\z')  from  (X2,  y2.  z2)  which  is 
the  location  of  the  "probe"  source,  and  z  =  Z2  *  z’.  Since  the  pressure  of  the  reversed 
wave  is  p3  =  Re[P3  exp(icot)],  inspection  of  (37)  shows  that  the  versed  wave  converges 
toward  S2  and,  if  it  is  allowed  to  pass  through  unobstructed,  diverges  when  it  passes 
beyond  S2.  Inspection  of  (36)  shows  that  the  wave  is  manifestly  a  conjugate  wave  where 
f  is  the  complex  "gain"  of  the  "phase  conjugate  mirror"  (PCM).  As  R  — >  0,  the  SPA 
breaks  down  and  it  is  necessary  to  include  the  effects  of  diffraction.  Let  d  denote  the 
smaller  of  the  diameters  of  the  bubble  layer,  the  pump  beam,  or  the  probe  beam.  By 
analogy  with  the  focal  properties  of  a  converging  lens  of  diameter  d,  it  may  be  argued  that 
the  maximum  pressure  at  the  focus  is  roughly 

(p3>m«  -  <d2/Xz2>  ^P2(Z  =  °)’-  (38) 

where  for  the  parabolic  approximations  used  in  the  derivations  of  Eq.  (32)  and  (33)  to  be 
applicable,  the  distance  to  the  focus  must  be  such  that  d  <  Z2.  It  is  remarkable  that  Eqs. 
(37)  and  (38)  give  the  same  expression  for  the  maximum  pressure  (apart  from  purely 
numerical  factors  close  to  unity)  as  the  appropriate  special  case  of  Eq.  (14)  of  Kustov  et 
al.36 

Sot  e  practical  consequences  of  this  analysis  are:  (i)  Since  PqC2  »  22,000  atm  for 
water,  generally  «  1  unless  lei  is  quite  large.  Large  values  of  lei  require  a  dense 
layer  of  bubbles  of  resonant  size.  For  a  dense  layer,  however,  it  becomes  necessary  to 
include  the  attenuation  of  the  pump,  probe,  and  reversed  wave  signals  which  was  not 
included  here,  (ii)  The  above  analysis  has  been  extended  to  allow  the  layer  to  be  at  a  plane 
Zq  shifted  from  the  plane  z  =  0.  This  introduces  a  phase  factor  exp(i  Iv^Zq/c)  in  P3  and 
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shows  that  while  adjacent  layers  of  bubbles  have  the  same  focal  point,  they  do  not 
necessarily  have  the  same  phase.  Hence  if  the  bubbles  in  a  layer  of  finite  thickness  L  are 
to  produce  reversed  waves  which  add  in  phase,  the  layer  should  be  sufficiently  thin  that 

co3L/c  «  k  which  gives  L  «  XJ2.  (39) 

For  a  thin  corrugated  layer,  L  becomes  the  maximum  amplitude  of  the  corrugations. 

This  calculation  is  significant  for  the  following  reason.  While  being  consistent  with 
results  in  Ref.  36,  the  analytical  methods  are  simpler  and  more  closely  parallel  than  those 
used  in  holography.  Equations  (37)-(38)  clarify  physically  relevant  magnitudes,  though 
we  make  use  of  previous  estimates  of  le(0)i,  o>2)l  for  bubbly  water  which  give38 

le(CD  ,0>2)l  -  nDX4,  D  *  10'2tol0‘\  (40) 

and  n  is  the  number  of  bubbles  per  volume  per  radius  size  increment  for  near  resonance 
sizes.  Inspection  of  Eqs.  (37)  and  (40)  shows  that  the  following  dimensionless  quantity 
will  need  to  be  large:  lelkL  =  nDX.3IV2lt.  The  analysis  leading  to  these  approximations 
was  needed  if  experiments,  described  below,  are  to  be  viable. 

B.  Experiments  on  Three-Wave  Mixing  in  a  Bubble  Layer 

A  graduate  student,  Steve  Kargl,  has  continued  to  study  the  generation  of  difference 
frequency  signals  using  a  modified  version  of  the  apparatus  shown  in  Fig.  20  of  the 
previous  Annual  Summary  Report .4  In  the  experiments  previously  described,  the 
nonlinear  layer  consisted  of  partially- stabilized  microbubbles  trapped  in  pores  (having 
diameters  of  a  few  pm)  in  a  sheet  of  Nuclepore  filter.  Scans  of  a  receiving  hydrophone 
revealed  localized  enhancements  of  a  difference  frequency  signal  which  were  suggestive  of 
the  focusing  of  a  reversed  wave.  The  frequencies  fj  =  C0j/27t  were  in  the  MHz  range. 
Subsequent  experiments  have  indicated,  however,  that  the  peaks  observed  were  probably 
not  due  to  a  reversed  wave.  Instead,  it  appears  that  the  linear  scattering  of  the  pump  and 
probe  waves  from  the  Nuclepore  sheet  (and  its  mount)  varies  with  hydrophone  position.  If 
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waves  at  frequencies  ©i  and  0)2  of  sufficient  intensity  are  incident  on  the  hydrophone,  a 
voltage  can  be  produced  at  the  difference  frequency  (03  =  0)1  -  a>2  as  a  consequence  of  the 
modulation  of  the  radiation  pressure  on  the  hydrophone;  this  can  mask  the  signal  of 
interest.  The  use  of  a  Nuclepore  filter  as  a  means  of  stabilizing  microbubbles  has  another 
drawback  which  we  were  not  aware  previously:  in  spite  of  the  hydrophobic  coating  on  the 
Nuclepore  filter,  the  gas  in  the  pores  can  dissolve  over  a  time  scale  of  a  few  hours  (or 
less).39 

Because  of  the  aforementioned  difficulties,  other  methods  have  been  (and  are  being) 
tried  for  creating  a  layer  of  bubbles.  The  bubbles  were  produced  by  electrolysis  and  lower 
frequencies  are  used:  typically  fi  *  400  kHz  and  f3  =  fi  -  f2  in  the  range  of  50  to  105 
kHz.  In  the  first  of  these  experiments  (spring,  1988)  the  electrolysis  wires  ran 
horizontally,  one  above  another,  in  an  array.  The  bubbles  came  off  in  roughly  vertical 
sheets.  The  sheets,  however,  were  neither  thin  nor  flat,  so  that  the  condition  given  as  Eq. 
(39)  is  typically  violated.  The  signal  amplitude  at  the  difference  frequency  varied  greatly 
from  burst  to  burst,  evidently  as  a  consequence  of  changes  in  shape  and  composition  of  the 
bubble  layer.  Hence,  while  interference  with  background  reflections  were  eliminated  with 
this  configuration,  the  time  variability  of  the  signal  ruled  out  the  use  of  simple  signal 
averaging  and  the  observations  were  inconclusive  as  to  the  focal  properties  of  the  reversed 
wave.  Instead  of  using  sophisticated  signal  processing  methods  to  look  for  evidence  of 
focusing  in  the  signal  at  f3,  an  alternative  arrangement  is  being  attempted  which  should 
facilitate  achieving  a  thin  stationary  layer  of  bubbles  as  well  as  a  low  background. 

Figure  14  is  a  diagram  of  the  approach  currently  being  tested.  Bubbles  are 
generated  by  electrolysis  at  a  horizontal  array  of  wires.  They  float  upward  and  are  trapped 
against  a  46  cm  x  46  cm  horizontal  sheet  of  clear  Mylar  having  a  thickness  of  12  pm.  The 
depth  of  the  500  gallon  tank  used  does  not  permit  the  sheet  to  be  in  the  far-field  of  the 
pump  source  at  S\  if  Si  is  directly  above  the  sheet.  Hence  the  beam  path  was  folded  by 
using  a  large  sheet  of  fused  silica  as  an  acoustic  mirror.  Taking  fi  =  400  kHz  and  f3  =  50 
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kHz  gives  focal  distances  from  Eq.  (35)  sufficiently  small  that  the  reversed  wave  comes  to 
a  focus  before  the  mirror.  The  pump  and  probe  amplitudes  at  the  bubble  layer  have 
magnitudes  of  103  Pa  and  102  Pa,  respectively  and  Eq.  (38)  suggests  that  the  IP3I  may  be 
as  large  as  1  Pa,  though  only  the  general  magnitude  of  the  bubble  density  has  been 
estimated.  Though  difference  frequency  signals  of  that  magnitude  have  been  observed,  it  is 
not  yet  known  if  they  are  focused.  An  effort  is  underway  to  automate  or  otherwise 
improve  the  speed  at  which  the  data  is  acquired.  This  is  necessary  since  there  are  changes 
in  the  bubble  layer  on  a  time  scale  of  a  few  hours  as  bubbles  dissolve  or  coalesce.  As 
expected  the  signal  at  f3  goes  away  after  several  hours,  evidendy  because  of  changes  in 
properties  of  the  bubble  layer. 

During  the  contract  period  under  review,  Kargl  demonstrated  a  simple  method  of 
locating  the  source  or  focus  of  a  diverging  or  converging  wave.  The  method  may  give  a 
definitive  way  of  demonstrating  the  focal  properties  of  a  reversed  wave,  though  at  the  time 
of  this  writing,  it  has  only  been  demonstrated  for  sound  received  direcdy  from  transducers. 
The  basic  idea  is  that  the  signal  is  recorded  for  various  hydrophone  positions  along  a 
straight  line.  Next  the  phase  is  determined  via  a  FFT  of  each  record.  From  the  phase  as  a 
function  of  position,  the  effective  source  location  is  determined. 

During  the  contract  period  under  review  the  experimental  facility  has  been  improved 
with  the  acquisition  of  a  power  amplifier  for  driving  the  probe  wave  source  and  a  digital 
delay  unit  for  timing  the  pump  and  probe  wave  bursts.  In  addition,  Kargl  has  calibrated  the 
principal  hydrophone  used. 

C.  Possibility  of  Phase  Conjugation  via  Reflection  from  a  Curved 
Vibrating  Surface 

In  a  previous  Annual  Summary  Report 21  the  possibility  of  phase  conjugation  via 
reflection  from  a  plane  vibrating  surface  was  discussed.  If  sound  of  frequency  {2  is 
incident  on  a  surface  which  vibrates  at  frequency  fi,  a  component  of  the  reflected  sound  is 


at  frequency  f3  =  f2  -  fi-  When  fi  =  2f2,  that  component  is  a  conjugate  wave  which 
converges  back  towards  the  source,  Fig.  15(a).  The  analysis  leading  to  that  result  uses  a 
paraxial  approximation  of  phase  shifts  and  hence  the  angle  of  incidence  relative  to  the 
surface  normal  is  not  large.  This  mechanism  was  independently  considered  in  the  Soviet 
literature.40 

During  the  contract  period  under  review,  Marston  reexamined  this  parametric 
mechanism  with  the  result  that  the  vibrating  surface  need  not  be  flat  For  example,  it  may 
be  a  section  of  a  sphere  or  a  cylinder  as  shown  in  Fig.  15(b).  The  essential  idea  is  that  each 
subsection  of  the  vibrator  could  be  treated  as  the  plane  reflector  in  Fig.  15(a).  The  use  of  a 
curved  reflector  has  at  least  two  distinct  advantages:  (i)  The  strong  signal  radiated  at 
frequency  fj  =  2f2  may  be  focused  away  from  the  spatial  region  of  interest,  and  hence 
reduce  the  background  sound  level  there,  (ii)  The  curved  reflector  may  be  made  part  of  a 
resonator  in  which  standing  waves  are  established.  Fig.  15(c).  Hence  the  energy  loss  may 
be  less  than  that  for  a  plane  open  vibrator  as  in  Fig.  15(a). 

Independently,  an  analysis  of  a  related  problem  recently  appeared  in  the  Soviet 
literature.41  Sound  of  frequency  f2  is  incident  (from  the  outside)  on  a  sphere  vibrating  at 
a  frequency  fj  =  2f2.  A  component  of  the  scattered  wave  is  predicted  to  be  a  conjugate 
wave  at  frequency  f2  directed  towards  the  source. 
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3.  Invited  address’,  P.  L.  Marston,  "Generalization  of  the  geometrical  theory  of  diffraction 
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conjugation  experiments:  Paraxial  analysis  and  analogy  with  optical  holography,"  [J. 
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